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We present simulations of atmosplieres of Earth-like aquaplanets that are tidally locked to their star, that is, planets whose orbital 
period is equal to the rotation period about their spin axis, so that one side always faces the star and the other side is always dark. 
As extreme cases illustrating the effects of slow and rapid rotation, we consider planets with rotation periods equal to one current 
Earth year and one current Earth day. The dynamics responsible for the surface climate (e.g., winds, temperature, precipitation) 
and the general circulation of the atmosphere are discussed in light of existing theories of atmospheric circulations. For example, as 
expected from the increasing importance of Coriolis accelerations relative to inertial accelerations as the rotation rate increases, the 
winds are approximately isotropic and divergent at leading order in the slowly rotating atmosphere but are predominantly zonal and 
rotational in the rapidly rotating atmosphere. Free-atmospheric horizontal temperature variations in the slowly rotating atmosphere 
are generally weaker than in the rapidly rotating atmosphere. Interestingly, the surface temperature on the night side of the planets 
does not fall below ~240 K in either the rapidly or slowly rotating atmosphere; that is, heat transport from the day side to the night 
side of the planets efficiently reduces temperature contrasts in either case. Rotational waves and eddies shape the distribution of 
winds, temperature, and precipitation in the rapidly rotating atmosphere; in the slowly rotating atmosphere, these distributions are 
controlled by simpler divergent circulations. The results are of interest in the study of tidally locked terrestrial exoplanets and as 
illustrations of how planetary rotation and the insolation distribution shape climate. 
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1. Introduction 

Planets generally evolve toward a state in which they be- 
come tidally locked to their star Torques the star exerts 
on tidal bulges on a planet lead to an exchange between 
the spin angular momentum of the planetary rotation 
and orbital angular momentum of the planet's revolution 
around the star, such that the rotation period around the 
spin a xis gradually ap proaches the orbital period of the 
planet (lHubbardlll984 . (The spin angular momentum of 
the star may also participate in this angular momentum 
exchange.) This process reaches its tidally locked end 
state when the rotation period is equal to the orbital pe- 
riod, so that one side of the planet always faces the star 
and the other side is always dark. The time it takes to 
reach this end state may exceed the lifetime of the plan- 
etary system, so it may never be reached (this is almost 
certainly the case for the Sun-Earth system). But plan- 
ets that are close to their star can reach a tidally locked 



state more quickly. Such close planets in other solar sys- 
tems are easier to detect than planets farther away from 
their star, and exoplanets that are believed to be tidally 
locked have indeed been detected in recent years (e.g., 
Charbonneau et al Here we investigate the atmo- 



spheric dynamics of Earth-like tidally locked aquaplan- 
ets through simulations with a three-dimensional general 
circulation model (GCM). Our purpose is pedagogic: we 
contrast rapidly and slowly rotating tidally locked Earth- 
like planets with each other and with Earth itself to illus- 
trate the extent to which atmospheric dynamics depend 
on the insolation distribution and planetary rotation rate. 

There are two areas of existing research on the at- 
mospheric dynamics of tidally locked planets. First, 
there are several studies motivated by "hot Jupiters" — 
large, clos e-in planets that have b een observed transiting 
a star (see S howman et al. I (l2009l) for a review). Second, 
there are studies of Earth-like exoplanets closely orbit- 
ing a relatively cool star — planets that have not yet been 
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but may be detected soon, for example, by NASA's re- 
cently launched Kepler space telescope or by ground- 
ba sed telescopes (as demon s trated by the obse rvations 
of ICharbonneau et al.ll2009l) . Ijoshi et all (1 19971) investi- 
gated the large-scale circulation of such Earth-like plan- 
ets and explored how th eir climate d epends on the mass 
of the atmosphere. And IJoshl ( 2003b documented their 
hydrological cycle using an Earth climate model. 

The studies by Joshi et al. provide a description of 
Earth-Uke tidally locked atmospheric circulations and 
how they depend on some parameters, such as the atmo- 
spheric mass. But some questions remain, among them: 
(i) How does the planet's rotation rate affect the circu- 
lation and climate? (ii) What controls the precipitation 
distribution (location of convergence zones)? (iii) What 
mechanisms generate large-scale features of the circula- 
tion, such as jets? (iv) What determines the atmospheric 
stratification? 

We address these questions by simulating tidally 
locked Earth-like aquaplanets with rotation periods equal 
to one current Earth year and one current Earth day. 
These are merely two illustrative cases: The slowly ro- 
tating case corresponds roughly to the tidally locked end 
state Earth would reach if the sun were not changing, 
that is, if the solar constant Sq = 1367 Wm~^ remained 
fixed. The rapidly rotating case corresponds to a terres- 
trial planet sufficiently close to a cool host star, so that 
the orbital period is one Earth day but the average inso- 
lation reaching the planet still is Sq = 1367 Wm^^, as 
presently on Earth. 

2. General circulation model 

We use an idealized atmosphere GCM with an active hy- 
drological cycle and an aquaplanet (slab ocean) lower 
boundary condition. The slab ocean has no explicit hor- 
izontal transports but is implicitly assumed to transport 
water from regions of net precipitation to regions of net 
evaporation, so that the local water supply is unlimited. 
T he GCM is a modified version of the model described 
in O' Gorman and Schneidei (l2008h. which is similar to 
the model described in Frierson et al. Briefly, the 

GCM is a three-dimensional primitive-equation model of 
an ideal-gas atmosphere. It uses the spectral transform 
method in the horizontal, with resolution T85, and fi- 
nite differences in cr = p/ps coordinates (pressure p and 
surface pressure Ps) in the vertical, with 30 unequally 
spaced levels. Subgrid-scale dissipation is represented 
by an exponential cutoff filter (iSmith et alji2002) . which 



acts on spherical wavenumbers greater than 40, with a 
damping timescale of 1 h on the smallest resolved scale. 
Most features of the simulated flows are similar at T42 
resolution or with different subgrid-scale filters; excep- 
tions are noted below. 

The GCM has a surface with uniform albedo (0.38) 
and uniform roughness lengths for momentum fluxes 
(5 X lO^'^m) and for latent and sensible heat fluxes 
(10^^ m). It has a gray radiation scheme, and a 
quasi-equilibrium moist convection scheme that relaxes 
convectively unstable atmospheric columns to a moist 
pseudo-adiabat with constant relative humidity 70% 
(lFrierson|[2007l) . Only the vapor-liquid phase transition 
is taken into account, and the latent heat of vaporization 
is taken to be constant. No liquid water is retained in the 
atmosphere, so any precipitation that forms immediately 
falls to the surface. Radiative effects of clouds are not 
taken into account, except insofar as the surface albedo 
and radiative parameters are chosen to mimic some of the 
global-mean effects of clouds on Earth's radiative bud- 
gets. 

Other model deta ils can be found in 
O' Gorman and Schneideii ( l2008h . However, 
the radiation scheme differs from that in 



O' Gorman and Schneideii (120081) and is described in 
what follows. 



a. Tidally locked insolation 

The top-of-atmosphere insolation is held fixe d at the in- 



stanta neous value for a spherical planet (e.g.. lHartmann 
1994 . 



'S'toa = 5*0 X max(0, cos(0) cos(A — Aq)) , (2.1) 

where (f> is latitude and A is longitude, with subsolar lon- 
gitude Ao = 270° and solar constant Sq = 1367 Wm"^. 
There is no diurnal cycle of insolation. That is, we as- 
sume zero eccentricity of the orbit and zero obliquity of 
the spin axis. 



b. Longwave optical depth 

O' Gorman and Schneideii (l2008l) used a gray radiation 



scheme with a longwave optical depth that was uniform 
in longitude and a prescribed function of latitude, thereby 
ignoring any longwave water vapor feedback. This is not 
adequate for tidally locked simulations with strong ra- 
diative variations in longitude. Therefore, we use a long- 
wave optical depth that varies with the local water vapor 
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concentration (cf. Thuburn and Crai jl2000l) . providing a 
crude representation of longwave water vapor feedback. 

As in Frierson et al. ( 2006), the longwave optical 
depth T has a term linear in pressure p, representing well- 
mixed absorbers like CO2, and a term quartic in pressure, 
representing water vapor, an absorber with a scale height 
that is one quarter of the pressure scale height. 



(2.2) 



Here, ti = 1.2 and po = 10^ Pa are constants. How- 
ever, the optical thickness of water vapor is a function of 
the model's instantaneous column water vapor concen- 
tration, 

1 fP" 

To = — qdp, (2.3) 
Pi Jo 

with specific humidity q and an empirical constant pi = 
980 Pa to keep tq order one for conditions typical of 
Earth's tropics. 

The details of the radiation scheme such as the con- 
stants chosen affect quantitative aspects of the simula- 
tions (e.g., the precise surface temperatures obtained), 
but not the large-scale dynamics on which we focus. The 
simulated su rface climate is similar to those presented in 
Joshil (12003 ), who used a more complete representation 
of radiation including radiative effects of clouds. This 
gives us confidence that the qualitative results presented 
here do not depend on the details of the radiation scheme. 

c. Simulations 

We conducted a rapidly rotating simulation with plan- 
etary rotation rate equal to that of present-day Earth, 
il = fls ^ 7.292 X 10~^s~\ and a slowly rotating sim- 
ulation with planetary rotation rate approximately equal 
to one present-day Earth year, = il£;/365. 

The results we present are averages over the last 
500 days of 4000-day simulations (with 1 day = 86400 s, 
irrespective of the planetary rotation rate). During 
the first 3000 days of the simulation, we adjusted the 
subgrid-scale dissipation parameters to the values stated 
above to ensure there was no energy build-up at the 
smallest resolved scales. 



L, is a measure of the importance of inertial accelera- 
tions (^C/^/i) relative to Coriolis accelerations {^fU) 
in the horizontal momentum equations. In rapidly ro- 
tating atmospheres, including Earth's, the Rossby num- 
ber in the extratropics is small, and the dominant balance 
is geostrophic, that is, between pressure gradient forces 
and Coriolis forces. In slowly rotating atmospheres, the 
Rossby number may not be small. If Ro = 0{l), inertial 
and Coriolis accelerations both are important, as in the 
deep tropics of Earth's atmosphere. If Ro 3> 1, Corio- 
lis accelerations and effects of planetary rotation become 
unimportant. In that case, there is no distinguished di- 
rection in the horizontal momentum equations, so the 
horizontal flow is expected to be isotropic, that is, the 
zonal velocity scale U and meridional velocity scale V 
are of the same order This is the dynamical regime of 
our slowly rotating simulation. 

In this dynamical regime, the magnitude of horizon- 
tal temperature variations can be estimated thro ugh scale 
analy sis of the horizontal momentum equations (ICharney 
19631) . In the free atmosphere, where frictional forces 
can be neglected, inertial accelerations scale advectively, 
like ~C/^/i, and are balanced by accelerations owing 
to pressure gradients, which scale like 6p/{pL), where 
Sp is the scale of horizontal pressure variations and p is 
density. The density and vertical pressure variations are 
related by the hydrostatic relation, p/H ^ pg, where 
H = RT/g is the pressure scale height (specific gas 
constant R and temperature T). If one combines the 
scalings from the horizontal momentum and hydrostatic 
equations, horizontal variations in pressure, density, and 
(potential) temperature (using the ideal gas law) scale 
hke 

Sp Sp 59 C/2 
p p e gH 



Fr, 



(3.1) 



where Fr = U^/{gH) is the Froude number For a 
terrestrial planet with H w 7 km and g = 9.8 ms^^, 
the Froude number is Fr < 10^^ for velocities U < 
25 ms^^. So free-atmospheric horizontal temperature 
and pressure variations are expected to be small insofar 
as velocities are not too strong (e.g., if their magnitude 
is limited by shear instabilities). This is the case in the 
tropics of Earth's atmosphere, and these expectations are 
also borne out in the slowly rotating simulation. 



3. Slowly rotating simulation 

The Rossby number Ro = U / (fL), with horizontal ve- 
locity scale U, Coriolis parameter /, and length scale 



a. Surface temperature 

In the slowly rotating simulation, the surface temperature 
mimics the insolation distribution on the day side of the 
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Figure 1 : Surface temperature (K) in simulations with e /365 (left) and fl e (right). The contour interval is 10 K. 



planet, decreasing monotonically and isotropically away 
from the subsolar point; the night side of the planet has 
a nearly uniform surface temperature (Fig. [T]i. Though 
the surface temperature resembles the insolation distri- 
bution, the influence of atmospheric dynamics is clearly 
evident in that the night side of the planet is considerably 
warmer than the radiative-convective equilibrium tem- 
perature of K. 

b. Hydrological cycle 

Surface evaporation rates likewise mimic the insolation 
distribution (Fig. |2]), as would be expected from a sur- 
face energy budget in which the dominant balance is 
between heating by shortwave radiation and cooling by 
evaporation. This is the dominant balance over oceans 
on Earth (e.g.. Trenberth et al. 2009) . an d in general on 
suffic i ently warm Earth-like aquaplane ts dPierrehumbert 
Gorman and Schneiden The dominance 

of evaporation in the surface energy budget in sufficiently 
warm climates can be understood by considering how the 
Bowen ratio, the ratio of sensible to latent surface fluxes, 
depends on temperature. For surface fluxes given by bulk 
aerodynamic formulas, the Bowen ratio, B, depends on 
the surface temperature, Tg, the near-surface air temper- 
ature, Ta, and the near-surface relative humidity, H, 



B 



L{q*{Ts)-nq*{Ta)y 



(3.2) 



where q* is the saturation specific humidity, Cp is the 
specific heat of air at constant pressure, and L is the la- 
tent heat of vaporization. Fig. [3] shows the Bowen ratio 
as a function of surface temperature, assuming a fixed 
surface-air temperature difference and fixed relative hu- 
midity. We have fixed these to values that are repre- 



sentative of the GCM simulations for simplicity, but the 
surface-air temperature difference and relative humid- 
ity are not fixed in the GCM. For surface temperatures 
greater than ~290 K, latent heat fluxes are a factor of 
> 4 larger than sensible h eat fluxes, as at Earth's sur- 
face jTrenberfli et al.ll2009h . Similar arguments apply for 
net longwave radiative fluxes at the surface, which be- 
come small as the longwave optical thickness of the at- 
mosphere and with it surface temperatures increase; see 
Pierrehumberj (l2009h for a more complete discussion of 
the surface energy budget. 

The precipitation rates also mimic the insolation dis- 
tribution (Fig.|2]l. There is a convergence zone with large 
precipitation rates (> 40 mm/day) around the subso- 
lar point. Precipitation rates exceed evaporation rates 
within ^-^15° of the subsolar point. Outside that region on 
the day side, evaporation rates exceed precipitation rates 
(Fig. |2]i, which would lead to the generation of deserts 
there if the surface water supply were limited. The atmo- 
spheric circulation that gives rise to the moisture trans- 
port toward the subsolar point is discussed next. 

c. General circulation of the atmosphere 

The winds are approximately isotropic and divergent at 
leading order They form a thermally direct overturning 
circulation, with lower-level convergence near the sub- 
solar point, upper-level divergence above it, and weaker 
flow in be t ween (Fig. |4]i. As in the simulations of 
Joshi et al.l (Il997h . the meridional surface flow crosses 
the poles. Moist convection in the vicinity of the sub- 
solar point results in strong mean ascent in the mid- 
troposphere there; elsewhere there is weak subsidence 
associated with radiative cooling (Fig. |5]). Consistent 
with the predominance of divergent and approximately 
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Figure 2: Evaporation (top row, mm/day), precipitation (middle row, mm/day), and evaporation minus precipi- 
tation (bottom row, mm/day) in simulations with fis/SGS (left) and He (right). Contours are shown at 0.1 x 
gO, 1,2,3,4 nimday"^ in the top two rows and —12.5, —1.0, 0, 1.0, 12.5 mm day in the bottom row. 




Figure 3: Bowen ratio vs surface temperature assuming a 3-K surface-air temperature difference, 70% relative humid- 
ity, and surface and surface air pressure of lO'^ Pa, using the same simpUfied saturation vapor pressure formula as in 
the GCM. 



isotropic flow, the Rossby number is large even in the extratropics (Ro > 10). 

The zonal mean zonal wind and streamfunction are 
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Figure 4: Zonal wind (left) and meridional wind (right) in il£;/365 simulation on the a = 0.28 (a), 0.54 (b), and 1.0 
(c) surfaces. The contour interval is 5 ms~^. 
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Figure 5: Pressure velocity (a; — Dp/ Dt) on the a = 0.54 surface in simulations with ri£;/365 (left) and il.E (right). 
Contours are shown at —0.2, 0, and 0.02 Pas^^. 



shown in Fig. |6] The zonal mean zonal wind is a weak residual of the opposing contributions from different lon- 
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Figure 6: Circulation in simulations with f7£;/365 (left panels) and He (right panels). Top row: zonal mean zonal 
wind (lines, contour interval 1 m ^ for ft e and 5 m ^ for Q,e/ 365) and divergence of the horizontal component of 
eddy (transient and stationary) angular momentum fluxes (colors, contour interval 1.0 x lO^^ms"^ for He/ 365 and 
1.5 X 10^^ ms^^ for He)- Bottom row: Eulerian-mean mass flux streamfunction (contour interval 25 x 10® kgs"^ 
for 17B/365 and 10 x 10® kg-s^^ for He). 



gitudes (Fig.|4]l. In the upper troposphere near the equa- 
tor, there are weak westerly (superrotating) zonal winds, 
as well as the eddy angular momentum flux convergence 
that, accord ing to Hide' s theorem, is necessary to sustain 
them (e.g., 'Hi delll96S " 
19801: [Schneidel|2006|) 



Schneider 1977: Held and Hou 



The Eulerian mean mass streamfunction consists of 
a Hadley cell in each hemisphere, which is a factor ^3 
stronger than Earth's and extends essentially to the pole. 
These Hadley cells are thermodynamically direct circu- 
lations: the poleward flow to higher latitudes has larger 
moist static energy {h = CpT + gz + Lq) than the near- 
surface return flow. The classic theory of Hadley cells 
with nearly inviscid, angular momentum-conserving up- 
per branches predicts that the Hadle y cell extent in- 



creases as the rotation rate decreases (Schneider 1977 



Held and Houlll980h . However, the results of this theory 
do not strictly apply here because several assumptions 
on which it is based are violated. For example, the sur- 
face wind is not weak relative to the upper-tropospheric 
winds (Fig. IDi, and the zonal wind is not in balance with 
meridional geopotential (or temperature) gradients (i.e., 



the meridional wind is not negligible in the meridional 
momentum equation). Nonetheless, it is to be expected 
that Hadley cells in slowly rotating atmospheres span 
hemispheres (IWilliamslll988l) . 



The instantaneous, upper-tropospheric zonal wind is 
shown in Fig.|7] and a corresponding animation is avail- 
able at http://www.gps.caltech.edU/~tapio/animations/tidallyJocked_slow.r 
Athough the flow statistics in the simulations are hemi- 
spherically symmetric (because the forcing and boundary 
conditions are), the instantaneous wind exhibits north- 
south asymmetries on large scales and ubiquitous vari- 
ability on smaller scales. The large-scale variability has 
long timescales (~80 days), and the zonal-mean zonal 
wind is sensitive to the length of the time averaging: for 
timescales as long as the 500 days over which we aver- 
aged, the averages still exhibit hemispheric asymmetries 
(Fig.©. 

The vertically integrated eddy kinetic energy is 8.4 x 
10^ J m"^ in the global mean. The instantaneous ve- 
locity fields exhibit substantial variability on relatively 
small scales (e.g.. Fig. |7]), and the kinetic energy spec- 
trum decays only weakly from spherical wavenumber 
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-55 -28 28 55 m s"^ 

Figure 7: Instantaneous zonal wind on cr 0.28 surface in the statistically steady state (on day 3900) of the two 
simulations. Animations are available at http://www.gps.caltech.edu/~tapio/animations/tidal lyjocked-slow.mov and 
|http://www.gps.caltech.edu/~tapio/animations/tidallyJocked_fast.mov| . 



~2G toward the roll-off near the grid scale owing to the 
subgrid-scale filter However, the peak in velocity vari- 
ance is at the largest scales (Fig.[8]l — as suggested by the 
hemispheric ally asymmetric large-scale variability seen 
in Fig.|7]and in the animation. 



d. Atmospheric stratification and energy transports 

Consistent with the scaling arguments presented above, 
horizontal variations in temperature are small in the 
free troposphere (above the ^750 hPa level. Fig. |9]l. 
This is reminiscent of Earth's tropics, where temperature 



variations are small because Coriolis accelerations are 



weak com^parej w i th ine rtial accelerations (e.g.,|Charney 
119631: LSobel et alJEx)!!) . Here, Coriolis accelerations 
are weak at all latitudes, and free-tropospheric horizontal 
temperature variations are small everywhere. 

The thermal stratification on the day side in the sim- 
ulation is moist adiabatic in the free troposphere within 
^30° of the subsolar point. Farther away from the in- 
tense moist convection around the subsolar point, includ- 
ing on the night side of the planet, there are low-level 
temperature i nversio ns (Fig. Wp), as in the simulations 
in [joshi et al. (1997). These regions have small optical 
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Figure 8: Spectrum of the mass-weighted, vertically averaged transient velocity variance for zonal (top row) and 
meridional (bottom row) wind components in the simulations with fl e /365 (left) and e (right). Contours are shown 
at 20,2\...,25m2s-2. 
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Figure 9: Temperature section of antisolar longitudes (a) and subsolar longitudes (b) of fi^ /365 simulation. Averages 
are taken over 10° of longitude. The contour interval is 10 K. 



thickness because of the low water vapor concentrations 
and therefore have strong radiative cooling to space from 
near the surface. The inversions arise because horizon- 
tal temperature gradients are constrained to be weak in 
the free troposphere, but near-surface air cools strongly 
radiatively, giving rise to inversions. 

The vertical and meridional integral of the horizontal 
moist static energy flux divergence as a function of lon- 
gitude is shown in Fig. [TO] (dashed curves). Energy di- 



verges on the day side of the planet and converges on the 
night side of the planet, reducing temperature contrasts. 
Near the subsolar point, there is substantial cancellation 
between the latent energy and dry static energy compo- 
nents of the moist static energy flux divergence as there 
is, for example, in Earth's Hadley circulation. As can be 
inferred from the minimal water vapor flux divergence on 
the night side (Fig. |2]l, the moist static energy flux diver- 
gence is dominated by the dry static energy (CpT + gz) 
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Figure 10: Left: Vertical and meridional integral of the divergence of the moist static energy flux, uh. Right: Dry 
static energy (red) and latent energy (blue) components of moist static enery flux divergence. Dashed lines for He/ 365 
simulation and solid lines for Q,e simulation. 



component on the night side (> 99% of the total). 



4. Rapidly rotating simulation 

In rapidly rotating atmospheres, the Rossby number in 
the extratropics is small, and geostrophic balance is the 
dominant balance in the horizontal momentum equa- 
tions. In the zonal mean, zonal pressure gradients van- 
ish, but meridional pressure gradients do not, so if the 
dominant momentum balance is geostrophic, winds are 
anisotropic U > V, where (•) denotes a zonal mean. 
Variations in the planetary vorticity with l atitude, /3, are 
central to v orticity mixing arguments (e.g., Rhine s|[l994 : 
|Held"2000^, which can account for the generation of at- 
mospheric jets: When a Rossby wave packet stirs the 
atmosphere in a region bounded by a polar cap, high- 
vorticity fluid moves equatorward and low-vorticity fluid 
moves poleward. This reduces the vorticity in the po- 
lar cap. By Stokes' Theorem, the reduced vorticity in 
the polar cap means the zonal wind at the latitude of the 
bounding cap decreases; if angular momentum is con- 
served, the zonal wind outside the polar cap increases. 
Irreversible vorticity mixing (wave-breaking or dissipa- 
tion) is necessary to maintain the angular momentum 
fluxes in the time mean. Thus, the larger planetary vor- 
ticity gradients of the rapidly rotating planet allow jets 
to form provided there is a source of wave activity that 
leads to vorticity stirring. 



We return to the analysis of lCharnevi ( 119631) to esti- 
mate temperature variations in the free atmosphere. For 
rapidly rotating planets, geostrophic balance holds in the 
horizontal momentum equations: 5p/{pL) ^ fU. Com- 
bining the scahng from the momentum equation with the 



hydrostatic relation, the pressure, density, and (potential) 
temperature variations scale like the ratio of the Froude 
number to the Rossby number. 



Sp 5p 
P P 



fUL 

gH 



Ft 
Ro' 



(4.1) 



Where the Rossby number is small (in the extratropics), 
the temperature variations will be a factor of order in- 
verse Rossby number (Ro^^ ^ 10) larger than in the 
slowly rotating simulation if similar values for U, g, and 
H are assumed. Thus, we expect larger horizontal tem- 
perature and pressure variations away from the equator 
in the rapidly rotating simulation. 

a. Surface temperature 

In the rapidly rotating simulation, the surface tempera- 
ture on the day side of the planet is maximal off of the 
equator and does not bear a close resemblance to the in- 
solation distribution; the night side of the planet has rel- 
atively warm regions in western high latitudes (Fig. [T]|. 
Compared to the slowly rotating simulation, the surface 
temperature is more substantially modified by the atmo- 
spheric circulation; however, the temperature contrasts 
between the day and night side are similar. 



b. Hydrological cycle 

Surface evaporation rates mimic the insolation distribu- 
tion (Fig. |2]i. This is one of the most similar fields be- 
tween the slowly and rapidly rotating simulations, as ex- 
pected from the gross similarity in surface temperature 
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and the smallness of the Bowen ratio at these tempera- 
tures (Fig. O. This might not be the case if the model 
included the radiative effects of clouds since the amount 
of shortwave radiation reaching the surface would be 
shaped by variations in cloud albedo which, in turn, de- 
pend on the atmospheric circulation. 

Precipitation rates are large in a crescent-shaped re- 
gion on the day side of the planet; the night side of the 
planet generally has small but nonzero precipitation rates 
(Fig. 121). The evaporation minus precipitation field has 
substantial structure: there are large amplitude changes 
from the convergence zones (P > E) to nearby areas of 
significant net drying {E > P). Comparing the slowly 
and rapidly rotating simulations shows that precipitation 
and E — P on the night side of the planet are sensitive to 
the atmospheric circulation. 

An interesting aspect of the climate is that the maxi- 
mum precipitation (near ^15° latitude) is not co-located 
with the off-equator maximum surface temperature (near 
^^40° latitude). The simulation provides an example 
of precipitation and deep convection that are not lo- 
cally thermodynamically controlled: the precipitation is 
not maximum where the surface temperature is max- 
imum; the column static stability (e.g.. Fig. \T3[ . and 
therefore the convective available potential energy, are 
not markedly different between the maxima in precipi- 
tation and temperature. However, if the surface climate 
is examined latitude-by-latitude instead of examining its 
global maxima, the region of large precipitation is close 
to the maximum surface temperature (as well as surface 
temperature curvature) at a given latitude. The structure 
of the surface winds, discussed next, and the associated 
moisture convergen ce are key for d etermining where pre- 



cipitation is large. ISobell (l2007h provides a review of 



these two classes of theories for tropical precipitation 
(thermodynamic control vs. momentum control) and the 
somewhat inconclusive evidence of which class of theory 
better accounts for Earth observations. 

The global precipitation is ~ 10% larger in the rapidly 
rotating simulation than in the slowly rotating simula- 
tion. This suggests that radiative-convective equilibrium 
cannot completely describe the strength of the hydrolog- 
ical cycle. 

c. General circulation of the atmosphere 

In the rapidly rotating simulation, the atmospheric circu- 
lation has several prominent features: there are westerly 
jets in high (^^65°) latitudes, the mid-tropospheric zonal 



wind exhibits superrotation, and the surface winds con- 
verge in a crescent-shaped region near the subsolar point 
(Fig.fTTTi. The equatorial superrotation and westerly jets 
are clear in the zonal mean zonal wind (Fig.|6]l. 

The existence of the high-latitude jets can be under- 
stood from the temperature field and eddy angular mo- 
mentum flux convergence (Figs. [13] and |6]l. There are 
large meridional temperature gradients, which give rise 
to zonal wind shear by thermal wind balance and pro- 
vide available potential energy for baroclinic eddies that 
transport angular momentum into the jets. In the vertical 
average, the eddy angular momentum transport into an 
atmospheric column is balanced by surface stress. Note 
that the region of higher temperatures on the night side of 
the planet coincides with the westerly jet's surface maxi- 
mum. The stronger winds at these latitudes lead to strong 
temperature advection. The equatorial superrotation is 
a conse quence of angular mor nentum flux convergence 
(Fig. leTl. lsaravanaiil jl993l) and lSuarez and Duffvl jl992h 
describe the emergence of superrotation generated by 
large-scale, zonally-asymmetric heating anomalies in the 
tropics. As in their idealized models, the zonal asymme- 
try in the low-latitude heating (in our simulation, pro- 
vided by insolation) generates a stationary Rossby wave. 
Consistent with a stationary wave source, in the rapidly 
rotating simulation, the horizontal eddy angular momen- 
tum flux convergence in low latitudes is dominated by the 
stationary eddy component. This aspect of the simulation 
is sensitive to horizontal resolution and the subgrid-scale 
filter. With higher resolution or weaker filtering, the su- 
perrotation generally extends higher into the troposphere 
and has a larger maximum value. 

There is a crescent-shaped region where the surface 
zonal wind is converging. This is where the precipita- 
tion (Fig. |2]i and upward vertical velocity (Fig. |5]l are 
largest. The horizontal scale of the convergence zone is 
similar to the equatorial Rossby radius, (c/ (3)^^^ ^ 10°, 
where /3 is the gradient of planetary vorticity and c is the 
gravity wave speed (estimated using a characteristic tro- 
pospheric value for the Brunt-Vaisala frequency on the 
day side of the planet). The surface zonal wind can be 
qualitatively understood as the equatorially-trapped wave 
response to stationary heating: equatorial Kelvin waves 
propagate to the east of the heating and generate easter- 
Ues; equatorial Rossby waves propagate to the west of 
the heating and generate westerlies (Gill 1980). 

The shape of the zero zonal win d Une and its hori- 
zontal scale are similar to those of the 'Gill' (mo) model, 
which describes the response of damped, linear shallow- 
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water waves to a prescribed heat/mass source. For the 
prescribed heat source in the original Gill model, the 
crescent-shape zero zonal wind line extends over ~2 
Rossby radii and, as in the GCM simulation, is displaced 
to the east of maximum heating on the equator 

A complicating factor in the analogy between the 
GCM's low-latitude surface winds and those of the Gill 
model is that the heating is prescribed in the Gill model, 
while it is interactive with the flow in the GCM. As previ- 
ously mentioned, the precipitation is strongly shaped by 
the winds. To see if the analogy between the winds in the 
GCM and in the Gill model breaks down because of the 
more complex structure of the latent heating, we force 
a variant of the Gill model with the GCM's precipita- 
tion field following the formulation of N eelin ( 1988 !) (see 
Appendix for details). The results of this calculation are 



compared with the GCM output in Fig. [12] The direction, 
large-scale structure, and, in the case of the zonal com- 
ponent, magnitude of the mass fluxes are similar between 
the GCM and precipitation-forced Gill model, though it 
is clear that there are quantitative differences. 

In contrast to the larger-magnitude zonal wind, the 
meridional wind is diverging at the surface and con- 
verging aloft near the subsolar point (right panel of 
Fig. [TTk.c). The dynamics of this are contained in the 
inviscid limit of the linear shallow-water equations: for a 
Sverdrup vorticity balance, the vortex stretching caused 
by the overall convergence near the surface near the 
equator must be balanced by moving poleward, toward 
higher planetary vorticity; hence the poleward merid- 
ional wind near the surface and the equatorward merid- 
ional wind aloft. 
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Figure 12: Zonal (left column) and meridional (right column) near-surface mass fluxes for the GCM (top row, aver- 
aged between the surface and a — 0.73 model level) and the Gill model forced by the GCM's precipitation (bottom 
row, see Appendix for Gill model equations and parameters). The contour interval is 2 x 10** kgm^^ s^^. 
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Figure 13: Temperature section of antisolar longitudes (a) and subsolar longitudes (b) of fi^ simulation. Averages are 
taken over 10° of longitude. The contour interval is 10 K. 



The Eulerian mean mass streamfunction (Fig. |6| has 
the opposite sense as Earth's Hadley cells: in the zonal 
mean, there is descent on the equator, poleward flow near 
the surface, ascent near 15°, and equatorward flow in the 
mid-troposphere. If the dominant balance in the zonal 
momentum equation is between Coriolis acceleration 
and eddy angular momentum flux divergence, —fv ^ 
—dyu'v' (i.e., small local Rossby number as defined in 



Walker and Schneider (l2006l) X then the angular momen- 
tum flux convergence that establishes the superrotating 
zonal wind also leads to equatorward mean meridional 
wind in the free troposphere, as in Earth's Ferrel cells 
(e.g.. lHddlEoO0l) . This can le ad to a dynamical fe edback 
that enhances superrotation dSchneider and Liul l2009i) : 
as superrotation emerges, the mean meridional circu- 
lation can change direction with a concomitant change 
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in the direction of mean-flow angular momentum fluxes 
(changing from exporting angular momentum from the 
deep tropics to importing it), which enhances the super- 
rotation. 

The instantaneous, upper tropospheric zonal 
wind in the rapidly rotating simulation is shown in 
Fig. I2I and a corresponding animation is available at 



Fig. [Tol l; there is substantial cancellation between dry 
static energy and latent energy flux divergence near the 
subsolar point. Though the hydrological cycle is more 
active on the night side of the rapidly rotating simulation 
than in the slowly rotating simulation, the dry static en- 
ergy component still dominates (> 80%) the moist static 
energy advection on the night side of the planet. 



|http://www.gps.caltech.edu/~tapio/animations/tidallyJocked_fastiM)^vo simulations are more similar in this respect 



The large-scale features of the general circulation such 
as the high-latitude jets and divergent zonal wind in the 
tropical upper troposphere are clear in the instantaneous 
winds. The eddies in the animation generally have larger 
spatial scales and longer timescales than in the corre- 
sponding slowly rotating animation. 

The vertically and globally integrated eddy kinetic 
energy is 1.0 x 10^ J m^^. This is about 20% larger 
than int the slowly rotating case. The eddy kinetic energy 
spectrum has a typical, steep shape to the small- 
est resolved wavenumber. In some sense, the integrated 
eddy kinetic energy belies the difference in synoptic vari- 
ability between the rapidly and slowly rotating simula- 
tions: in the extratropics, for zonal wavenumbers be- 
tween ^3-6, the rapidly rotating simulation has a fac- 
tor of 2-3 times more velocity variance (Fig.[8]l than the 
slowly rotating simulation. 

d. Atmospheric stratification and energy transports 

The tropospheric temperatures on the day side of the 
planet (Fig. \13[ reflect the surface temperature varia- 
tions: the temperature field has a local maximum near 
'--^40° and is relatively uniform up to high latitudes (pole- 
ward of ^50°). In the free troposphere on the day side, 
the lapse rates are close to the moist adiabatic lapse rate, 
computed using the local temperature and pressure, over 
a region roughly within the 300 K contour of the surface 
temperature. Note that 300 K does not have a particular 
physical significance — we are simply using it to describe 
a feature of the simulation. There is a local minimum of 
temperature on the equator which may be related to the 
downward vertical velocity there (Fig. ID. In low lati- 
tudes, there is a near-surface inversion on the night side 
of the rapidly rotating simulation that, as in the slowly ro- 
tating simulation, is the result of weak temperature gra- 
dients in the free troposphere and the small optical thick- 
ness of the atmosphere. 

As in the slowly rotating simulation, the moist static 
energy flux divergence is positive on the day side and 
negative on night side of the planet (solid curves in 



than might have been anticipated given the differences 
in their flow characteristics, although there are regional 
differences that are obscured by averaging over latitude. 
The broad similarities can be understood by considering 
the moist static energy budget. In the time mean, denoted 
(•), neglecting kinetic energy divergence and diffusive 
processes within the atmosphere, the mass-weighted ver- 
tical integral, denoted (•), of the moist static energy flux 
divergence is balanced by surface energy fluxes, Fg, and 
radiative tendencies, Qiad^ 

V-(T[7^)=F3 + 0,,d. (4.2) 

As a result of the gross similarity of the two simulations 
in evaporation and low-latitude stratification (in part due 
to the smaller dynamical role that rotation plays near the 
equator) and hence the radiative cooling, the divergence 
of the moist static energy flux is also similar In the 
extratropics, considering the sources and sinks of moist 
static energy does not provide a useful constraint because 
the stratification is dy namically determined (e.g., iHeld 



1978t ISchneiderll2007l) and there can be geostrophically 



balanced temperature gradients in rapidly rotating atmo- 
spheres; as result, the radiative cooling, through these dy- 
namical influences on temperature, is determined by the 
flow, so the moist static energy flux divergence will not in 
general be independent of rotation rate. Indeed, the warm 
regions on the night side of the rapidly rotating simula- 
tion have larger moist static energy flux divergence than 
in those regions on the night side of the slowly rotating 
simulation as a result of the difference in circulation. 

5. Conclusions 

We have examined the dynamics of tidally locked. Earth- 
like planets. In particular, the dynamical regime of the 
atmosphere depends on the planet rotation rate, as antic- 
ipated from scaling arguments. The simulations demon- 
strate the importance of the atmospheric circulation in 
determining the surface climate. For example, the differ- 
ence between the precipitation distributions in the slowly 
and rapidly rotating simulations clearly shows the depen- 
dence of the hydrological cycle on the circulation of the 
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atmosphere. Interestingly, some aspects of the simulated 
climate are not sensitive to the planet's rotation rate. In 
particular, the temperature contrast between the day and 
night side and evaporation rates are similar between the 
slowly and rapidly rotating simulations. 

The general circulation of the slowly rotating atmo- 
sphere features global-scale, thermodynamically-direct, 
divergenct circulations with ascending motion at the sub- 
solar point and descending motion elsewhere. In con- 
trast, the general circulation of the rapidly rotating atmo- 
sphere features extratropical jets that owe their existence 
to rotational (Rossby) waves, and tropical surface winds 
that are the result of stationary equatorial (Rossby and 
Kelvin) waves. The isotropy of the winds in the slowly 
rotating regime and the anisotropy of the winds in the 
rapidly rotating regime are expected from the dominant 
balance of the horizontal momentum equations. Expec- 
tations for the free-atmospheric pressure and temperature 
gradients based on scale analysis were realized in the 
simulations of the different circulations regimes: tem- 
perature gradients are weak where Coriolis accelerations 
are weak (in low latitudes of the rapidly rotating atmo- 
sphere and globally in the slowly rotating atmosphere) 
and larger where Coriolis accelerations are dynamically 
important (in the extratropics of the rapidly rotating at- 
mosphere). 

While many aspects of the simulations can be ex- 
plained by general circulation theories for Earth and 
Earth-like atmospheres, there remain quantitative ques- 
tions that require more systematic experimentation than 
we have attempted here. For example, the dependence of 
the surface climate on the solar constant may be differ- 
ent from simulations with Earth-like insolation because 
of the differences in stratification. 
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6. Appendix: Gill model forced by precipitation 

We use th e formulation for a Gill model forced by pre- 
cipitation ("N eelinl 1988 1). as opposed to the standard mass 
sink. The model equations are 

-fv = -d^cf) - e^u (6.1) 

fu = ~dy(i> - £mV (6.2) 

(?{dxU + dyv) = —aP — £t<P, (6.3) 

where a = LR/cp is the combination of the latent heat 
of vaporization, gas constant, and heat capacity at con- 
stant pressure that results from converting the precipi- 
tation rate to a temperature tendency in the thermody- 
namic equation and approximating Newtonian cooling 
by geopotential damping, e„i is the Rayleigh drag co- 
efficient of the momentum equations, is the damp- 
ing coefficient for the geopotential, other variables have 
their standard meanings (e.g. P is precipitation), and 
(•) — J^^° (•) dp/g is the mass-weighted vertical inte- 
gral over the boundary layer (pb is the pressure at the top 
of the boundary layer). We use fairly standard values for 
the constants (Table [T]i. The qualitative aspects of the so- 
lutions that we discuss are not sensitive to these choices. 

The three equations can be combined into a single 
equation for v: 

[-eriel + /?'y')/c' + + dyy) + P = 

-{e-nidy - Py d-x) aP/c^, 
(6.4) 

where /3-plane geometry (/ = (3y) has been assumed. 
The boundary value problem for v is discretized by 
Fourier tranforming in the x-direction and finite differ- 
encing in the y-direction. It is solved using the GCM's 
climatalogical precipitation field on the right-hand side 
(Fig.13. 
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constant 


value 


c 


50ms-i 


P 


2.22 xlO-" 




(1 day)-i 




(20 day)-i 


a 


7.14 xlO^ Jkg-i 



Table 1: Gill model constants and their values. 



Frierson, D. M. W., 2007: The dynamics of idealized convec- 
tion schemes and their effect on the zonally averaged tropical 
circulation. / Atmos. Sci, 64, 1959-1976. 

Frierson, D. M. W., I. M. Held, and P. Zurita-Gotor, 2006: A 
gray-radiation aquaplanet moist GCM. Part I: Static stability 
and eddy scale. /. Atmos. Sci., 63, 2548-2566. 

Gill, A. E., 1980: Some simple solutions for heat induced trop- 
ical circulation. Quart. J. Roy. Meteor. Sac, 106, 'XAl-A&l. 

Hartmann, D., 1994: Global Physical Climatology. Academic 
Press, 411 pp. 

Held, I. M., 1978: The tropospheric lapse rate and climatic sen- 
sitivity: Experiments with a two-level atmospheric model. J. 
Atmos. Sci., 35, 2083-2098. 

Held, I. M., 2000: The general circulation of the atmosphere. 
Proc. Program in Geophysical Fluid Dynamics, Woods 
Hole Oceanographic Institution, Woods Hole, MA. URL 
[https ://darchive.mblwhoilibrary.org/handle/ 1912/15 

Held, I. M., and A. Y. Hou, 1980: Nonlinear axially symmetric 
circulations in a nearly inviscid atmosphere. J. Atmos. Sci., 
37, 515-533. 

Hide, R., 1969: Dynamics of the atmospheres of the major 
planets with an appendix on the viscous boundary layer at 
the rigid bounding surface of an electrically-conducting ro- 
tating fluid in the presence of a magnetic field. / Atmos. Sci., 
26, 841-853. 

Hubbard, W. B., 1984: Planetary Interiors. Van Nostrand Rein- 
hold Co., 334 pp. 

Joshi, M., 2003: Climate model studies of synchronously rotat- 
ing planets. Astrobiology, 3, 415-427. 

Joshi, M. M., R. M. Haberle, and R. T. Reynolds, 1997: Simu- 
lations of the atmospheres of synchronously rotating terres- 
trial planets orbiting M Dwarfs: Conditions for atmospheric 
collapse and the implications for habitability. Icarus, 129, 
450-465. 



Neelin, J. D., 1988: A simple model for surface stress and low- 
level flow in the tropical atmosphere driven by prescribed 
heating. Quart. J. Roy. Meteor Soc, 114, 1^1 -HQ. 

O'Gorman, P. A., and T. Schneider, 2008: The hydrological cy- 
cle over a wide range of climates simulated with an idealized 
GCM. / Climate, 21, 3815-3832. 

Pierrehumbert, R. T, 2002: The hydrologic cycle in deep-time 
climate problems. Nature, 419, 191-198. 

Pierrehumbert, R. T, 2009: Principles of Planetary Climate. 
Cambridge University Press, 531 pp. 

Rhines, P B., 1994: Jets. Chaos, 4, 313-339. 

Saravanan, R., 1993: Equatorial superrotation and maintenance 
of the general circulation in two-level models. / Atmos. Sci., 
50, 1211-1227. 

Schneider, E. K., 1977: Axially symmetric steady-state models 
of the basic state for instability and climate studies. Part II. 
Nonlinear calculations. / Atmos. Sci., 34, 280-296. 

Schneider, T, 2006: The general circulation of the atmosphere. 
Ann. Rev. Earth Planet. Sci., 34, 655-688. 

Schneider, T, 2007: The thermal stratification of the extratrop- 
ical troposphere. T. Schneider and A. H. Sobel, Eds., The 
Global Circulation of the Atmosphere, Princeton University 
Press, 47-77. 

Schneider, T, and J. J. Liu, 2009: Formation of jets and equa- 
torial superrotation on Jupiter. / Atmos. Sci., 66, 579-601. 

Showman, A. P, J. Y.-K. Cho, and K. Menou, 2009: Atmo- 
spheric circulation of extrasolar planets. S. Seager, Ed., Ex- 
oplanets. University of Arizona Press. 

Smith, K. S., G. Boccaletti, C. C. Henning, I. N. Marinov, C. Y. 
Tam, I. M. Held, and G. K. Vallis, 2002: Turbulent diffusion 
in the geostrophic inverse energy cascade. / Fluid Mech., 
469, 13-48. 

Sobel, A. H., 2007: Simple models of ensemble-averaged trop- 
ical precipitation and surface wind, given the sea surface 
temperature. T. Schneider and A. H. Sobel, Eds., The Global 



JAMES-D 



Tidally locked Earth 



17 



Circulation of the Atmosphere, chap. 8, Princeton University 
Press, 219-251. 

Sobel, A. H., J. Nilsson, and L. M. Polvani, 2001: The weak 
temperature gradient approximation and balanced tropical 
moisture waves. J. Atmos. Set., 58, 3650-3665. 

Suarez, M. J., and D. G. Duffy, 1992: Terrestrial superrotation: 
A bifurcation of the general circulation. / Atmos. Sci., 49, 
1541-1554. 

Thubum, J., and G. C. Craig, 2000: Stratospheric influence on 



tropopause height: The radiative constraint. J. Atmos. Sci., 
57, 17-28. 

Trenberth, K. E., J. T. FasuUo, and J. Kiehl, 2009: Earth's 
global energy budget. Bull. Amer. Meteor. Soc, 90, 31 1-323. 

Walker, C. C., and T. Schneider, 2006: Eddy influences on 
Hadley circulations: Simulations with an ideaUzed GCM. 

J. Atmos. Sci., 63, 3333-3350. 

Williams, G. P., 1988: The dynamical range of global circula- 
tions — I. Climate Dyn., 2, 205-260. 



Journal of Advances in Modeling Earth Systems - Discussion 



